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Abstract. - We propose a simple configuration of a one-dimensional parallel array of Josephson- 
junctions in which the pinning potential for trapped fluxons lacks inversion symmetry (ratchet 
potential) . This system can be modelised by a set of non-linear pendula with alternating lenghts 
and asymmetric harmonic couplings. We show, by molecular dynamics simulation, that fluxons 
behave as single particles in which the predictions for overdamped thermal ratchets can be easily 
verified. 



Directional motion (DM) of brownian particles has been a target of basic and applied 
research in the last five years. The initial motivation and interest in this field came from cell 
biology: The study of the mechanism of vesicles transport inside eukariotic cells, via motor 
proteins along microtubules 0, ^. Later on new systems with the same underlying ideas for 
transportation were proposed. Those systems include: phase separation engines||^, drop mo- 
tion under ac forces growth of surfaces ||^, and rectification in asymmetric superconducting 
rings (SQUID) To have DM of a particle submitted to a periodic one-dimensional potential 
V{x) = V{x + L) it is needed i) broken spatial symmetry i.e. V{x) ^ V^—x) and ii) to drive 
the particle out of thermal equilibrium. This last condition makes possible to extract work 
from the system. 

The first condition is usually fulfilled by the use of an "ad hoc" ratchet potential [0 . The 
latter condition can be implemented in different ways: by an oscillating driving field , by a 
time-correlated non-thermal noise Q or a fluctuating V{x) potential Q, among others. In all 
these cases, the conditions for the fluctuation-dissipation theorem do not hold. 

In this paper we propose a new experimental realization of DM in a very simple and 
controllable system: a parallel Josephson-junction array. Such system has become an excellent 
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Here is the figure 1 



Fig. 1. - Schematic equivalent circuit of the parallel array of Josephson- junction. External current 
is injected from top to bottom. Rs is a small shunt resistance to deal with overdamped junctions. 



realization of the Frenkel-Kontorova (FK) or discrete sine-Gordon model Fluxons (or 
kinks) in the array move along a periodic pinning potential (the so called Peierls-Nabarro 
potential) . We will show that in a simple geometry this potential is a ratchet one and a fluxon 
behaves as a particle with well defined directional motion. 

This Letter is organized as follows: we begin by describing the basic assumptions for the 
array and the equations of motion for the phases are obtained. Then, we calculate the energy 
profile for a fiuxon trapped in the array and show its "ratchet" character. Finally, we compare 
our results for the dynamics of the fluxon with what is expected from single particle dynamics. 

We shall study first the array configuration. Let us consider a parallel array of Josephson- 
junctions with alternating critical currents {Id and Ic2, generally Ica) and alternating area 
plaquettes. The area differences cause alternating self-inductions Li and In order to deal 
with overdamped junctions, we shunt each of them with a small external resistance Rg so that 
we can neglect capacitance effects. Figure 1 is a schematic picture of the circuit. We denote 
by 2'!TXj the phase difference across the junction j. I{t) is an external driving current. The 
equation for current conservation through the junction j read as 

^ + Ica sin 27rxj = /(i) + - /* = I{t) + l]_^ - l] (1) 

Jig 

where 

V, = ^oxj (2) 
being $o is the fiux quantum. The phase diference along a plaquette is given by 

- = 1^ (3) 

where 

i„ (/*+/]') (4) 

is the magnetic fiux through a plaquette, sum of two contributions: external and induced 
magnetic fields]^. We set now 

^ext ^ 0. Combining equations (|l|), (|), (|) and (§, and 
dividing the array in two sublattices, the equations of motion for the phases are: 

+ Id sin 2Trxj ^ I{t) + ^i^j-i - ^j) + ^(^J+i ~ ^j) 
— + Ic2 sin27ra;j+i = I{t) + ^{xj - xj+i) + ^(^3+2 - Xj+i) (6) 
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Normalizing the equations by II — ^o/2Li and denoting hy K = 47r/ciii/$o, ct = Ic2/Ici 
and /3 = L1/L2 we obtain the dimensionless equations of motion 

K 

Xj + — sin27ra;,- — I{t) + (a^j-i — Xj) + (3{xj^i — xj) (7) 
Xj+1 + sin27ra;j+i = i{t) + I3{xj - Xj+i) + {xj+2 - xj+i) (8) 

ZTT 

where time derivatives are given in the time scale r — 2Li/Rs. Here K play the role of 
discreteness parameter Q and a and (3 measure the asymmetry of the array. Flux trapped 
in the array determines the boundary conditions, therefore, if we have n flux quanta in N 
plaquettes, then xn+i = Xi+n. Note that these equations of motion (|^) and (||) correspond 
to the well known FK model in a generalised asymmetric form. The overdamped dynamics for 
the symmetric model have been reviewed in reference . The mechanical analogue of the J J 
array we are dealing with is a set of pendula with two different lengths coupled by two kinds 
of harmonic springs. 

It is expected for if — > that the pinning potential for fluxons vanishes. But we are 
interested in the opposite limit in which fluxons are strongly pinned to the discrete lattice. 
Due to the spatial asymmetry of the array we expect to have an asymmetric pinning potential. 

Now we will try to describe the fluxon as a single particle in the pinning periodic potential. 
Fluxon location is given by its center of mass (CM) which can be expressed as ||l^ 

1 ^ 

Xcm{xi, ■ ■ ■,xn) = - + ^i{xi+i - x^) (9) 

1=1 

for any given phase configuration {xi}. On the other hand, the energy of the array is the sum 
of two contributions: Josephson and magnetic energy. In the adimensional parameters defined 
above this energy reads as 



E{xi,...xn)^ ^ 



— — (cos27ra;i+ a cos 27r.Ti+i) + -(xi+i - a; J + - [x,+2 - x.^ 
47r"^ 2 2 



i—l,odd 



(10) 

From these expressions we can define E{Xcm) — rnin^^.yE{xi . . . xn) such that {xi} are kink 
configurations whose CM is Xcm- 

In order to calculate this potential, several approaches have been tried in the context of 
non- linear discrete lattices. For low pinning (low K), a good estimation is obtained using the 
collective coordinates method |l2]. Essentially this method assumes a soliton profile for the 
{xi} corresponding to the continuous limit (sine-Gordon) and, in some cases, an analytical 
solution can be reached which only converges for low values of K. We have adopted a more 
general way to (numerically) obtain the potential profile E{Xcm)- Maxima of this energy 
correspond to saddle-points in the N-dimensional phase space. For configurations containing 
one fluxon, N-1 directions are stable and one is unstable. Such maxima points can be obtained 
using standard minimax methods. We have performed a linear stability analysis around these 
points to get the direction of destabilization of the maximum energy configuration. Using 
this saddle configuration as initial condition, we perturb it along the unstable direction and 
study the relaxation (by numerical integration of equations (|^) and (|^), with I{t) — ) to the 
adjacent minima. Then, the energy and CM along the trajectory are computed (according to 
equations (^ and 10). Using this method we have been able to calculate the energy profile 
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Here is the figure 2 



Fig. 2. - Energy profile E{Xcm) for different values of asymmetry parameters {a, 13) (a) a = 0.5, {3 — 
1.0. Same plaquette areas and different critical currents, (b) a = l,/3 = 0.5. Different plaquettes and 
the same critical currents, (c) a = 0.5, f3 = 0.5 gives "ratchetlike" potential. 

E{Xcm) for a trapped fluxon as a function of its CM. Figure 2 shows potential profiles for a set 
of the model parameters {K, a, (3). In the following we set K = 4.0. For a = 0.5 and (3 = 1.0 
(different critical currents), energy profile shows a double well structure symmetric respect to 
the bottom of the wells. The a = 1.0, (3 = 0.5 (different areas) potential is symmetric respect 
to the tops. As expected, for any other values the potential profiles do not show inversion 
symmetry. The values a = 0.5 and P — 0.5 give a good approximation to the asymmetric 
sawtooth potential used in the literature It is interesting to note that, unlike other works 
in extented systems p^ , p^ , neither the on-site potential nor the interparticle potential are 
asymmetric in the field variables. 

We will study now the dynamical behavior of a fluxon in the asymmetric lattice. For all 
simulations, we take N = 30, n = 1, K = 4.0, a = 0.5 and /3 — 0.5. We drive the system out of 
thermal equilibrium by applying an external ac bias currents. The i(t) term in equations (|^) 
and (^ is then expressed as 

i{t) = iacSinLOt + ^{t). (11) 

Here, ^{t) is white noise {{^{t)) = and {^{t)C{t')) = 2Dd{t - i'))0 which, in absence of 
other forces, brings the system to thermal equilibrium. Thus, the equations of motion take 
the form of a system of stochastic differential equations. We have solved them using a fourth 
order Runge-Kutta method for the deterministic part and third order for the stochastic one 

[0. 

We will concentrate first on the deterministic (T — 0) dynamics. We check the ratchet 
behavior measuring the critical current of the array under dc driving. For positive driving we 
find a fluxon deppining current /J^ w 0.18 whereas for negative « 0.31 (for the symmetric 
case, a — 1.0 and /3 — 1.0, Idp ~ 0.19). Under ac currents we should observe a IV curve 
showing rectiflcation of the external current and the maximum efficienty for rectification (in 
the low frequency limit) of the ac currents is expected to be around /j^ 1^. Figure 3 shows 
the IV (voltage versus amplitude lac) curves for Idc = and T = for different frequencies 

(^) Where D = Jlf^r the temperature in units of the magnetic energy. 
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Here is the figure 3 



Fig. 3. - Voltage-Current (IV) curves at T = and ac current, for different frequencies u> in units of 
1/r. (a) Lu = 27r0.01 . (h) cu ^ 27r0.05. (c) uj = 27r0.1. 



oj. Low frequency response clearly resembles to that found for a single particle ||T^. As it was 
noted in one particle simulations, voltage is quantised: 

V= — 2^{x;) = i-uj), (12) 

being p and q integer numbers [ p6| . For finer resolution, the IV curves appears to have a devil's 
staircase structure (see inset of figure 3c) . This quantisation result can be straightforward 
obtained by assuming that, in the overdamped ac dynamics, the phases can be expressed 
by a two dimensional envelope function |l^ Xi{t) = ({i — Vt/^o,ujt) [Q. Now, imposing 
energy balance along the trajectory 

V / x^idt^Y] / Xii{t)dt (13) 
Jo Jo 

we find that only voltages commensurate with the current frequency give non-zero contribution 
to the right side of eq. [l^. 

Quantised steps appears rounded at non-zero temperature as a consequence of multiples 
hops between dynamically close attractors, corresponding to the mode-locking trajectories |1^. 
In the adiabatic (low frequency) and low temperature limit we have fitted the IV curves with 
the solution of Magnasco |^ for a sawtooth ratchet potential. We have found a good fit (see 
figure ^ for asymmetry ratio A1/A2 = 0.56 and energy barrier AE — 0.16 obtained from 
the computed energy profile (see figure 2c). At low frequencies the steps disappear, even for 
very low temperatures. At high frequency the steps seem to be more stable against thermal 
fluctuations and probably could be observed in real experiments. 

In summary, we have proposed a new experiment for DM based in the motion of a fluxon in 
a one-dimensional Josephson-junction array. We have calculated the effective pinning potential 
for the fluxon and found the parameters needed to have almost a sawtooth one. The simulated 
IV curves mimic the behavior of a single particle in an asymmetric periodic potential. 

The actual feasibility to fabricate these kinds of arrays along with the simplicity of our 
design make this experiment really accesiblc pO|. Evenmore, the possibility of introducing in 
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Here is the figure 4 



Fig. 4. - Voltage-Current (IV) curves for ac current at D = 0.001 and uj — 27r0.01. Thick line shows 
the comparison with equation 3 of reference Q. 

the array a controlled number of fluxons could serve to study the influence of interaction in 
directional motion brownian particles Preliminary results in this directions show that 

complex behavior dominates the multifluxon dynamics. 
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